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Abstract 

Mechanosensation is a key part of the sensory repertoire of a vast array of different 
cells and organisms. The molecular dissection of the origins of mechanosensation 
is rapidly advancing as a result of both structural and functional studies. One 
intriguing mode of mechanosensation results from tension in the membrane of the 
cell (or vesicle) of interest. The aim of this review is to catalogue recent work that 
uses a mix of continuum and statistical mechanics to explore the role of the lipid 
bilayer in the function of mechanosensitive channels that respond to membrane 
tension. The role of bilayer deformation will be explored in the context of the 
well known mechanosensitive channel MscL. Additionally, we make suggestions 
for bridging gaps between our current theoretical understanding and common 
experimental techniques. 
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1 Mechanosensation and the Channels that Me- 
diate it 

Cells interact with each other and with their external environment. These inter- 
actions are enabled by transmembrane proteins - machines that have evolved to 
allow cells to detect and respond to changes in their environment. These proteins 
detect external cues, such as an increase in ligand concentration or the presence of 
forces or voltage, and transiently alter the permeability of the cell membrane allow- 
ing ions, water, or even larger mo l ecules to cross as wel l as triggering receptors for 
signaling (iBarry and Lynch! 120051 : IClapham et all . 1200 ll ) . The passage of these ions 
(or molecules) and the triggering of receptors then leads to a series of downstream 
events within the cell, enabling a response to these environmental cues. 

Mechanical forces and their corresponding deformations constitute one of the 
most important classes of external cues. Mechanosensation is a widespread phe- 



nomenon in a host of dif f erent single-celled and multicellular organisms (jFain . 



2001 



20031: iGillespie and Walkerl l200ll; Ikatsumi et al 12004 iKloda and Martinac 
Nauli and Zhou . 2004 ; Sachsl . 1991 ). In bacteria, experimental evidence suggests 
that mechanosensation arises to detect and regulate the response to changes in the 
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osmotic environment (jChang et all . Il998l ; iPivetti et all . 120031 ; ISukharev et all . Il997l ). 
More generally, the issue of cell shape and its attendant deformation is impor- 
tant not only in the context of osmotic str ess and the management o f physical 
stresses to which membranes are subjected (IMorris and Homannl . l200ll ). but also 
arises i n context of remodeling of the cell and organelle me mbranes during cell di- 
vision ( jChristensen and Strangd . 1200 ll ; iKamada et all . Il995l ) . 

In multicellular organisms, mechanosensation is important in a variety of ways. 
One intriguing class of mechanosensors is linked to motility. For example, in ne- 
matodes like the much studied C. elegans, mechanosensation permits the worm to 
decide which way to move and may have a role in de tecting body curvature, thu s 



telling the worm when to change its wave-like shape (I Gillespie and Walker . l200ll ). 
Similarly, flies have hair bristles that respond to touch (IDuggan et all . 120001 ) . while 
the mechanosensitive lateral-line organelles in zebrafish provide the means for de- 
tecting di rectional water movement in a way very similar to the workings of our 
inner ear (iGillespie and Walked . 1200 ll ). In each of these cases, genetics has led to 
the identification of a variety of genes implicated in the ability of the organism to 
respond to some form of mechanical stimulus. Parallel insights have been obtained 
in plants (Arabidopsis in particular), with the ide ntification of a collection of noy el 
proteins that also appear to be mechanosensitive (IHaswell and Meyerowitzl . 120061 ). 

Mechanosensitive ion channels are a class of membrane proteins that have re- 
cently garnered significant interest. Genetic, biochemical and structural studies all 
conspire to make this a particularly opportune time to demand a more quantitative 
picture of the function of these channels. In particular, there is a growing list of 
success stories in which the structures of channels associated with mechanosensation 



have been found in both closed and op e n sta tes (IBass et all . 120021 ; IChang et all . 11998 



Perozo et all . l2002al : IPerozo and Reed . 120031 ) . In addition, functional studies that 



probe how gating depends upon membrane tension or external forces are beginning 
to mak e it possible to dissect the various contributions to the energetics of channel 
gatin g (jAkitake et all . 120051 ; IChiang et all . |2004 IPerozo et all . l2002bl ; ISukharev et all . 



19991 ). 



As a result of these studies, a number of ideas have been proposed to explain the 
different ways in which external force can couple to membrane-protein conformation. 
Two modes of action that have been hypothesized for channels are: i) cases in which 
physical, polypeptide linkers pull on some part of the protein resulting in gating, ii) 
cases in which tension in the surrounding bilayer forces the channel to open. The 
aim of this article is to show how statistical mechanics and simple models of bilayer 
elasticity can be used to glean insights into this second class of mechanosensors. 

The remainder of the article is built in four main sections. In the next section, 
we describe how statistical mechanics can be used to analyze the probability 
that a two-state mechanosensitive channel is open. This discussion will include 
an analysis of how the external load (i.e. the tension) can be included in the 
statistical mechanical treatment of these problems. The next section considers the 
elastic deformations imposed on a bilayer by the presence of a transmembrane 
protein, and shows how these deformations result in a mechanosensitive channel 
acting as a bistable switch (i.e. a protein with two stable conformations). In the 



2 



THE ROLE OF LIPID BILAYER MECHANICS IN MECHANOSENSATION 



subsequent section, we discuss experimental considerations that will help form 
a tighter connection between theory and experimental techniques. Finally, we 
examine the way multiple channels in a membrane might interact through the 
intervening lipid bilayer and how these interactions can alter the conformational 
statistics of individual channels. 



2 Statistical Mechanics of Mechanosensitive 
Channels 

To begin, we review the application of statistical mechanics to a simple two-state 
mechanosensitive channel. This analysis will serve as the starting point for our sub- 
sequent, more detailed analysis which explores how bilayer elasticity can contribute 
to the energetics of the closed and open states of a channel. 

2.1 Lipid Bilayer vs. Protein Internal Degrees of Freedom 

One convenient scheme for characterizing the state of ion channels is to invoke the 
state variable a, which is defined by o = if the channel is closed and a = 1 if 
the channel is open. Our aim is to compute the open probability P open which, in 
terms of our state variable a, can be written as (a), where (• • • ) denotes an average. 
When (a) ~ 0, this means that the probability of finding the channel open is low. 
Similarly, when (a) « 1, this means that it is almost certain that we will find the 
channel open. To evaluate these probabilities we need to invoke the Boltzmann 
distribution, which tells us that the probability of finding the system in a state 
with energy E(a) is p(a) = e~@ E ^ / Z , where Z is the partition function defined by 
Z = e~^ E ^\ (3 = l/ksT, ks is Boltzmann's constant, and T is the temperature 
in degrees Kelvin. 

On the level of a single channel, we introduce e c iosed and e open for the energies 
of the closed and open states, respectively, as shown in Figure [H These energies 
contain contributions from deformations of the surrounding lipid bilayer as well 
as internal protein energetics; however, they do not contain the tension-dependent 
driving force which we will address separately. The state variables can be used to 
write the channel energy (in the absence of tension) as 

E{a) = (1 - <j)e dosed + ae open . (1) 

With these energies in hand, we can assign weights to the different states as shown 
schematically in Figure [TJ Within this scheme, the probability that the channel is 
open is given by (cr) and can be computed as (a) = J2t=o a P( a )-> where p(a) is the 
probability of finding the channel in state a. To compute these probabilities, we 
invoke the Boltzmann distribution, and evaluate the partition function given by 

i 

cr=0 
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STATE ENERGY WEIGHT PROBABILITY 




Figure 1: States, Boltzmann weights and corresponding probabilities for a two-state 
ion channel. The two different states have different energies and the probability of 
these different states is determined by the Boltzmann distribution. 



As a result, we see that the open probability can be written as 

^ ^ g P e closed -(- (D — 0£open ^ _J_ f^P( e open~^closed) ^ ^ 

This expression is relatively sterile in the absence of some term that tunes the en- 
ergies of the open and closed states to reflect the impact of external driving forces. 
In fact, one of the most remarkable features of ion channels is that the probability 
of being in different states can be tuned by external factors such as ligand concen- 
tration, the application of a voltage, or application of tension in the surrounding 
membrane. In general, this formalism can account for any of these driving forces, 
but we will restrict our attention to the important case of mechanosensitive channels, 
where the key driving forces are mechanical. In this case, gating occurs when the 
energy balance between the open and closed states is altered by membrane tension. 

To give the origin of membrane tension a physical meaning, we introduce the 
notion of a "loading device" , which we define as the external agent acting on a lipid 
bilayer to alter its tension. As depicted by hanging weights on the bilayer in Figure 
[2J we can make a toy model of how changes in bilayer geometry are coupled to the 
energy of this loading device. The point of introducing this hypothetical situation is 
to enforce the idea that, in our statistical mechanical treatment of this problem, the 
loading device is an important part of the overall free energy budget of the system. 
As a result, when we write down the partition function for a problem involving a 
channel and a deformable membrane, we have to account for the internal protein en- 
ergetics, the deformation energy of the lipid bilayer, and the energy associated with 
the loading device itself. In particular, we note that an increase in the membrane 
area will lead to a lowering of the weights depicted in Figure [2] and a corresponding 
decrease in the energy of the loading device. Of course, the application of tension in 
real membranes is not perf ormed by hanging weights, but through techniques such 
as micropipette aspiration ( Goulian et al . 1998 ; Rawicz et al . 2000l ). Nevertheless, 



the concept of hanging weights brings the importance of the energy of the loading 
device into sharp focus. 
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Figure 2: Energy of the loading device for membrane deformation. This figure com- 
pares the unloaded and loaded membrane and shows how membrane deformation 
results in a lowering of the potential energy of the loading device. In this hypothet- 
ical experiment, the tension (force per unit edge length) in the membrane is given 
by r = mg/Al where Al is the distance between two consecutive hooks, and g is 
the acceleration due to gravity. 



For the case of tension-activated ion channels, the open probability, (a), is dic- 
tated by a competition between the energetic advantage associated with reduction 
in the energy of the loading device and the energetic cost of the open state due to 
both the internal protein energetics and the energetics of membrane deformation. 
Following up on the idea of Figure [2j but now with special reference to the case 
of a mechanosensitive ion channel, Figure [3] shows how the opening of the channel 
results in a reduction of the energy of the loading device. 

The total area of the bilayer is constant (to within a few percent), and as a result, 
when the channel opens and the radius gets larger the weights in our hypothetical 
loading device are lowered by some amount, which lowers the potential energy. The 
greater the weights, the larger the change in potential energy. The notion of weights 
is a simple representation of externally applied forces on the membrane. If we 
imagine a finite membrane with fixed £1X6 cl cLS shown in Figure [3J, when the channel 
opens, the outer radius will change as AR out = (R/ R out )AR, where R is the closed 
channel radius, AR is the change in channel radius upon opening, R out is the outer 
radius of the membrane when the channel is closed, and AR out is the increase in 
the outer radius of the membrane when the channel opens. We are interested in 
evaluating the change in potential energy of the loading device {i.e. the dropping 
of the weights) as a result of channel opening. To do so, we compute the work 
associated with the force F, which is most conveniently parameterized through a 
force per unit length (the tension, r) acting through the distance AR out as shown 
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Figure 3: Schematic of how channel opening results in a relaxation in the loading 
device. For simplicity, we represent the loading device as a set of weights attached 
to the membrane far from the channel. When the channel opens, these weights are 
lowered, and the potential energy of the loading device is decreased. 



in Figure [3J This results in 



R . ~ S^Rout 



AG tension = tAs X — AR X - . (4) 

. tiout v 

force on arc s v y ^— v— 

displacement of arc number of arcs 

where AG represents a change in free energy. We have introduced the variable As 
for the increment of arc length such that r = F/ As. Given these definitions, we see 
that the change in the energy of the loading device is given by 

AGtension = ~t2tcRAR. (5) 

In light of our insights into the energy of the loading device, we introduce the 
energy as a function of the applied tension r, which is given by 

E(a, r) = (1 - a)e dosed + ae open - otAA. (6) 

The term — otAA favors the open state and reflects the fact that the energy of 
the loading device is lowered in the open state. In fact, this term reveals that any 
increase in protein area is energetically favored when membrane tension is present, 
which could imply hidden mechanosensitivity in other classes of ion channels and 
receptors - a subject we will touch upon later in this review. 

To compute the open probability of the channel in the presence of applied tension, 
we need to once again evaluate the partition function Z = ^2 a e~P E ( a \ Using the 
energy given in eqn. [6j we find 

g — ^-fit-closed _|_ e -P(e open -TAA) ^y-j 

This permits us to write down the open probability directly as 

„-f3(e pen-TAA) 1 
p = C = (g\ 

open e -f3(e open -rAA) _|_ e -/3e closed ]_ _|_ e /3(e p e „-e ciosed -rAA) ' V > 
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STATE ENERGY WEIGHT PROBABILITY 




Figure 4: States, weights and corresponding probabilities for a two-state 
mechanosensitive channel under load. 



The corresponding states, weights, and probabilities for a channel under applied 
tension are shown in Figure HI The open probability of a mechanosensitive channel 
is shown in Figure as an increasing function of the applied tension. 

To understand how a particular channel is going to behave under a driving force, 
we need to know two things. First, we need to understand the channel's intrinsic 
preference for each of its two states, which is encoded by e c i OS ed and e open . Second, 
we need to understand how the external driving force alters the relative energies of 
these different states. With these two quantitative measurements in hand, statistical 
mechanics allows us to compute the behavior of the channel under a range of driving 
forces. To make further progress, we need to examine the microscopic origins of tdosed 
and e open . Intriguing recent experiments suggest that these energies are driven in 
large measure by membrane deformations. 



3 Bilayer Free Energy and Gating of a 
Mechanosensitive Channel 

The abstract formalism of the previous section leaves us poised to examine 
mechanosensation to the extent that we can understand the physical origins of 
^closed and e open . The main idea of this part of the review is to show how sim- 
ple models of the elastic properties of lipid bilayers can be used to determine the 
bilayer's contribution to e c \ ose d and e open . One of the key clues that hints at the 
importance of membrane deformation in dictating channel gating is the data shown 
in Figure In particular, this plot shows how the open probability depends upon 
the lipid carbon tail length. This data strongly suggests that the energetics of the 
surrounding membrane is an important part of th e overall free energy budget of 



channel gating (also see iMartinac and Hamilll (120021 ) and the informative review by 
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Figure 5: Ion channel open probability as a function of applied tension. The plot 
shows P op en = (°") as a function of the applied tension r. The parameters used 
in the plot for a model mechanosensitive channel are e open — e c i ose d = 10 k^T and 
AA = 10 nm 2 . The critical tension is LOA^T/nm 2 , corresponding to P ope n = 1/2 . 
For reference, the tension can be rewritten as 1 pN/nm ~ 0.25 /c^T/mri 2 . 



Jensen and Mouritsen j200i ). 



The parameters e c iosed and e open can each depend on some combination of the 
energetics of protein conformation, membrane deformation, and hydration energy. 
Our strategy is to use the tools of continuum mechanics to calculate how the defor- 
mation of lipids surrounding a protein and the app lied tension work in concert to 
affect the channel's preference for a particular state ( Dan and Safranl . Il998l ; iHuang , 



19861 ; iNielsen et all . Il998l ; IWiggins and Phillips! . 120051 ). Unfortunately, relatively lit- 



tle is known about how the internal rearrangements of the protein and the hy- 
dration energy of the channel pore contribu te to the overall free energy balance 
(lAnishkin et all . 120051 ; lYoshimura et all . |l999j). This ignorance is in part due to a 
lack of general rules that tell us how internal rearrangements translate into changes 
in protein energy. Further, the lack of crystal structures in the open and closed 
states of many channel proteins means we cannot be sure where each residue moves, 
which are exposed to the surrounding lipids and which are facing the hydrated in- 
ternal pore. It is also difficult for molecular dynamic s to comment on the energies 



associated with the internal movements of the protein (lElmore and Dougherty! . 120011 . 
20031 : IGullingsrud et all . l200ll ; IGullingsrud and Schultenl . l2Q03h because the all-atom 



energies of these simulations are very large in comparison to the changes in free en- 
ergy, and hence it is difficult to distill relatively small free energy changes in the 
background of large energy fluctuations. To complicate the issue further, it is also 
possible that the internal movements of the protein yield relatively small free energy 
changes between the two states, but may provide various kinetic hurdles in the form 
of energy barriers, which affect the transition rate from one state to another. 

It is reasonable on the scale of a single membrane protein to ask whether a 
bilayer composed of discrete lipid molecules can be approximated as a continuum 
material. We argue heuristically that, gi ven the relat i ve di f fusion coefficients of 
membrane proteins (D ~ 0.1 — lfim 2 /s) (IDoeven et all . 120051 ; iGambin et all . 120061 ; 
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Figure 6: Ion channel open probability for different lipids. The graph shows the 
open probability of the mechanosensitive channel MscL as a function of the applied 
pressure across the bilayer for three different lipid tail lengths. Pressure difference 
is related to bilayer tension via a constant, and hence this suggests that bilayer 
thickness (for carbon tail l e ngths of 16, 18 and 20) affects channel function. Figure 
adapted from lPerozo et all (j2002bl ). 



Guigas and Weisd . 120061 ) and lipids (D ~ 10/im 2 /s) (IKahya et all . 120031 ) , in the time 



it takes a transmembrane protein to diffuse one lipid diameter, many lipids will have 
exchanged places near the protein, in a sense averaging out the discreteness of the 
lipid mol ecules. Additionally, the trans ition time for protein conformational change 
(~ 5 (is) (jShapovalov and Lesterl . 120041 ) is slow compared to lipid diffusion. Hence, 
we argue the bilayer can be approxi mated as a continuou s material in equilibrium 
with well-defined elastic properties (IHarroun et all . |l999j). Further, we choose to 
formulate our an alysis in the la nguage of continuum mechanics, rather than lateral 
pressure profiles (jCantorl . Il999n . 

Ap proximating t h e membrane a s a continuum mate ri al dDan and Safran j. 



1998; 



Harrou n et all. Il999l; iHelfrichl . Il973l ; iHuangj . Il986l ; iNielsen et al 



1998 



Wiggins and Phillips! . 120041 . 120051 ) , we will concentrate our analysis on how the me- 



chanical properties and deformations of lipids affect the energy balance of the pro- 
tein, and how tension can play the role of a driving force for gating the channel. In 
particular, the mechanosensitive channel of large conductance (MscL) is one of the 
best characterized mechanosensitive channels. Additionally, a combination of X-ray 
crystallography and electron paramagnetic resonance studies h aye yielded insights 



into the structures of bot h the closed and open states of MscL (IChang et all . 11998 



Perozo et all . l2002al . l200ll ). One of the outcomes of this structural analysis is the idea 



that the structure can be roughly approximated as a cylinder, making it amenable 
to mechanical modeling. MscL exemplifies many of t he characteristics one might 
call " design principles" for a mechanosensitive channel (jWiggins and Phillips! . 12004 . 



20051 ). such as change in hydrophobic thickness, a change in radius, and sensitivity 



to membrane curvature. In the remainder of the review, we will lay the foundation 
for a continuum mechanical understanding of how lipid deformations and tension 
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work together to give a switchable channel. 



3.1 The Case Study of MscL 



In the prokaryotic setting, the physiological purpose of MscL i s thought to be an 



emergency r e lief yalve under cond i tions of hypoosmotic shock (jChang et all . 11998 



Pivetti et all . 120031 ; ISukharev et all . Il997l ). whereby the osmotic pressure difference 



between the inside of a cell and the environment translates into increased mem- 
brane tension. The channel responds by gating and non-selectively releasing os- 
molytes to the environment u ntil the internal and external pressures are equilibrated 



(ISukharev et all . 12001 
First 



19991 ). This presents us with (at least) two key questions, 
what gives MscL its ability to "sense" tension in the membrane? Second, 



what role is the lipid bilayer playing in the gating transition? 

We will argue that the answers to these questions are found in the properties of 
a lipid bilayer and the geometrical features of the channel as revealed in Table HJ 
In particular, the bilayer has four key elastic properties that give it the ability to 
transduce tension and resist deformation by a transmembrane protein. The most 
striking elastic feature is the in-plane fluidity of the bilayer, which, in the absence 
of cytoskeletal interactions, results in equalization of tension throughout the mem- 
brane. This means that any in-plane stress (i.e. tension) on the membrane is felt 
everywhere equally. Hence, in the case of MscL, an increase in tension is applied 
uniformly to the outer edge of the protein, essentially trying to "pull" the channel 
open. We argue it is this "pulling" which constitutes the driving force for channel 
gating. However, this driving force is competing with the energetic cost to gate the 
channel due to internal conformational changes within the protein and deformations 
of the surrounding lipid. 

Three other properties give the membrane the ability to store energy elasti- 
cally upon deformation. First, each leaflet of the membrane resists changes in the 
angle b etween adjacent lipid molecules, lead i ng to bending s tiffness of the mem- 



brane (|Dan and Safran . 1998 ; lHarroun et all. 19991 ; iHelfrichl . 1 19731 ; iHuangj . |1986 



Nielsen et all . Il998l ; I Wiggins and Phillips! . 120051 ). Second, the membrane has a pre- 
ferred spacing of the lipid mole cules in-plane and will resist any change s in this 
spacing due to external tension (IDan and Safran! . 1 19981 ; iRawicz et all . 120001 ). Third, 
the membrane has a well-defined equilibrium hydrophobic thickness which, when 
given an embedded protein of a differ ent hydrophobic thickness, leads to energet 



1998; Harroun et al. 1999 



ically costly ' hydrophobic mismatch' (IDan and Safran 
Nielsen et al Il998l : Iwiggins and Phillips! . I2004J . l2005h . 

The competition between the driving force and the energetic cost to gate the 
channel hints at a set of design principles that dictate how the channel behaves as 
a bistable switch. If we neglect the molecular details of MscL, its conformational 
change can be characterized by a set of simple changes in geometrical parameters. 
In particular, in our coarse-grained description we will think of the gating transition 
as being accompanied by changes in height, radius and protein angle, all of which 
couple to various modes of membrane deformation as shown in Figure [71 The central 
question becomes, is deformation of the lipids surrounding the protein a major 
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rdl dlilctcl . 


Value; 


oouice. 


OlObeQ iieigllt 


o.o rim 


onang et ai ^ryyoj 


Closed radius 


2.5 nm 


oUang et al (lyyoj 


Open height 


2.5 nm 


rcrozo et ai ^zuuzaj 


Open radius 


3.5 nm 


rcrozo et ai (zuuzaj 


Measured ZA/i 


zu nm 


omang et ai ^zuu4j 


Measured AG* 


51A; B T 


Chiang et al (2004) 


Calculated AG* (at critical tension) 




this article 


Critical Tension* 


~ z.o K^jJ /nm 


omang et ai (zuu4j 


ijytic leiision 


~ o.o rugi i nm 


nawicz et ai ^zuuuj 


Bending Modulus (Kb) 


~ 20 A; B T 


Niggemann et al (1995) 
Rawicz et al (2000) 


Area Stretch Modulus (K A ) 


~ 60 k B T/nm 2 


Rawicz et al (2000) 


Leaflet Thickness (7) 


1.75 nm 


Rawicz et al (2000) 



Table 1: MscL geometrical and bilayer elastic parameters. (*> These parameters 
depend on the elastic properties of the bilayer, in particular the bilayer bending 
modulus (k & ), the bilayer area stretch modulus (Ka), and the leaflet hydrophobic 
thickness (/). 



player in gating energetics? Indeed, experiments have already suggested that the 
gating characteristics are intimately linked to the hy drophobic mismatch b e tween 
the protein and bilayer as was shown i n Figur e [6] (IJensen and Mouritsenl . 12004 ; 
Martinac and Hamuli . 120021 ; iPerozo et all . l2002bl ). It is the goal of the following 



sections to build up a theoretical framework for understanding the various kinds 
of bilayer deformation around a transmembrane protein and to describe how these 
deformations contribute to the overall free energy budget associated with the gating 
of MscL (and probably other channels as well). 



3.2 Bilayer Deformation, Free Energy and the Role of Ten- 
sion 

To investigate the contribution of membrane deformation to channel gating in 
mechanosensitive channels, we put our ignorance of the internal protein energet- 
ics aside and focus on the response of the membrane. The point of this analysis 
is to see how large the membrane contributions are to the free energy of channel 
gating, and to examine how they compare to the measured values. A mechanosen- 
sitive channel must resist the driving force due to tension to exhibit the properties 
of a bistable switch. As we will demonstrate in this section, deformation of the 
surrounding lipids can provide this resistance, and almost certainly does in the case 
of MscL, given our knowledge of the open and closed structures and th e body of ex- 
perimental data describ i ng the interactions b etween lipids and MscL (IPerozo et all . 
2002bl ; IPowl et al 120031 : ISukharev et~al Il999t ). 
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The deformations that a transmembrane protein induces can be most broadly 
split into two main classes: those that deform the midplane of the bilayer, and those 
that deform the bilayer leaflet thickness. If the deformation i s not too severe, these 
two types of deformation are independent of one another (j Wiggins and Phillips! . 
20051 ). Figure [7] shows these two classes of deformation and the simple model ide- 
alizations implied by elastic descriptions. The basic structure of the models we 
consider are those in which the contributions of deformation to the overall free en- 
ergy are obtained by computing local bending and thickness deformation, and then 
summing over the contributions from all the area elements making up the bilayer. 



3.2.1 Midplane Deformation 

Deformation of the midplane of the bilaye r involves a c ost to bend th e mid 



plane from its flat, equil ibrium position (lHarroun et all . Il999l ; iHelfrichl . Il973l ; 
Wiggins and Phillips! . 120051 ). We use the function h(r) to denote this change in 
height of the bilayer midplane as a function of the position r as shown in Figure [71 
The energy cost associated with bending the membrane away from its flat configu- 
ration can be written as 



G 



(mid) 
bend 



Kb 

1 



J (V 2 Mr)-C Q ) 2 d 2 r, 



(9) 



where the b ilayer bending modulus Kb — 20 ksT ( iNiggemann et all . 11995 
Rawicz et all . |2000| ) and C Q is the midplane spontaneous curvature. Throughout 



the review the gradient operator is defined by V = (d/dx, d/dy), and the Laplacian 
operator by V 2 = d 2 /dx 2 + d 2 /dy 2 , in Cartesian coordinates. In general, bilay- 
ers with symmetric leaflet compositions have zero midplane spontaneous curvature. 
Tension also plays a role in the energetics of midplane deformation because any bend 
in the midplane results in a reduction in the projected area of the membrane, which 
couples directly to an increase in the energy of the loading device. This effect is 
quite intuitive when one considers deformations of a macroscopic membrane under 
tension and results in a contribution to the free energy of the form 



G ( mid ) = 1 

ten 2 



(V/i(r)) 2 d 2 r 



(10) 



where the tension, r, ranges from zer o up to the nominal membrane lytic tension of 
~ 3.5 /c^jT'/nm 2 Q flRawicz et all . l2000h . In general, the elastic parameters we use are 
representative of a typical phosphatidylcholine (PC) lipid. Thus the total energy 
expended to deform the midplane over an area A is 



G ( mid ) 



^g(VMr)) 2 + |(V 2 Mr)) 2 ) d 2 r. 



(11) 



The logic behind this kind of analysis is to find the free energy minimizing function 
h(r). One way to carry out this minimization is by solving a partial differential 



*The lytic tension of a bilayer is technically a dynamic quantity ( Evans et all . 12003). however, 
we quote the lytic tension as the tension at which bilayer lysis is a rapid, spontaneous process. 
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equation that is generated by formally minimizing the free energy. An alternative 
(and approximate) scheme to be explored later in this section is to make a guess 
for the functional form of h(r) and to minimize with respect to some small set of 
parameters. This approach is called a variational method and can be quite useful 
for developing intuition. 

In the midplane-deforming model, the protein can dictate the slope of the mem- 
brane at the protein-lipid interface which, in addition to the protein radius, will 
determine the deformation energy. The length scale over which the membrane re- 
turns to its unperturbed state is given by \fi%jr and the energy for this type of 
deformation is 



G (mid) (i?,r) = 6 2 irR y 



Kb 



(12) 



where R is the radius of the protein, 9 is the slope of the membrane at the protein- 
lipid interface as show n in Figure [71 and Kj are modified Bessel function s of the 
second kind of order i ([Turner and Send . 120041 ; IWiggins and Phillips! . 120051 ) . Given 
a protein with a particular radius and fixed boundary slope, an increase in tension 
will make any deformation more costly. Hence, for midplane deformation, increased 
tension prefers a flatter membrane and/or smaller protein radius. To get a feel for 
the energy scale of this deformation several examples for different parameter values 
are summarized in Table [2j 

With the contribution to the free energy difference arising from midplane defor- 
mation in hand, we can now explore the competition between applied tension and 
the energetics of membrane deformation in dictating channel gating. The key to 
understanding the interplay between tension and deformation energetics lies in the 
scaling of these two effects with protein radius. The midplane deformation energy 
scales roughly linearly with the radius of the protein and is unfavorable. On the 
other hand, the term proportional to the applied tension scales as the square of the 
protein radius and favors the open state. If we fix the membrane slope, then the 
energy of a midplane deforming protein as a function of protein radius and tension 
is 

G(R,t)^G^{R,t)- rnR? . (13) 
membrane loading device 

As tension increases, the potential energy of the loading device will eventually over- 
come the deformation energy and a larger protein radius will be the preferred state. 
Indeed, midplane deformati ons have been hypothe sized to be an important func- 
tional mechanism of MscL ( Turner and Send . 12004 ) . One of the uncertainties that 
accompanies a model of this type is the fact that there is some function that connects 
the slope of the membrane at the protein-lipid interface (8) with the current radius 
of th e channel, that is, there is some unknown function 9(R) ([Spencer and Reed . 
2002j). Future experiments will be necessary to further clarify this point. If we 
make the simplest approximati on that 6(R) = constant and look at two reasonable 
values of 9 = 0.6 and 6 = 0.8 (ITurner and Send . 12004 ). using eqn. [13] and the pa- 
rameters in TableHJ we find the rather small critical tensions ~ 0.004 /c^T/mri 2 and 
~ 0.06 /c^T/mri 2 , respectively, compared to the known critical tension of MscL at 
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Figure 7: Modes of bilayer deformation. Differences in the equilibrium shape of 
the membrane and an embedded protein give rise to local deformations, a) The 
undeformed state is a transmembrane protein with zero hydrophobic mismatch and 
a flat midplane. Even an initially undeformed membrane exhibits tension depen- 
dence since tension induces bilayer thinning, b) An angled protein induces midplane 
bending, characterized by the function h(v) and the boundary slope 9. As tension 
increases, the most preferred energetic state is 9 = 0. c) A membrane protein that is 
thinner than the equilibrium thickness of the bilayer compresses the bilayer causing 
local area expansion and bending of each leaflet, characterized by the function u(r). 
d) A membrane protein that is thicker than the equilibrium thickness of the bilayer 
stretches the bilayer causing local area reduction and bending of each leaflet, e) An 
increase in tension will decrease the energetic cost of a membrane protein that is 
thinner than the equilibrium thickness of the bilayer, as the membrane thins and 
approaches zero hydrophobic mismatch, f) Likewise, an increase in tension will in- 
crease the energetic cost of a membrane protein that is thicker than the equilibrium 
thickness of the bilayer. 
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Tdble 2: Typicdl free energies for midplane deformation. The first row indicdtes 
how tension leads to an increase in deformation energy. The second and third rows 
show the sensitivity to the boundary slope. The fourth row indicates how protein 
radius changes deformation energy. The last row is a comparison with the known 
radius change and critical tension of MscL. 



~ 2.5 k B T/nm 2 ( Anishkin et all . 2005 ; Chiang et al . 2004 ). Though we have shown 
that midplane deformations are capable of endowing a channel protein with bista- 
bility, the scale of the critical tension and the free energy difference between confor- 
mations indicates that, at least for MscL, an additional kind of deformation might 
be important as well. 



3.2.2 Thickness Deformations 



We have examined how protein conformation can alter midplane bending of the 
surrounding lipid bilayer and how this deformation energy penalizes the open state 
by virtue of its larger radius. A second major class of deformations are those that 



bend and compress a single leaflet of the membrane (Aranda-Espinoza et all . |1996 



20051 ) and can be thought of 



Huang . 1986 : Nielsen et all . 1998 : Wiggins and Phillips , 
as imposing a local thickness on the lipid bilayer that is different from its equilibrium 
value, as illustrated in Figure [71 This kind of deformation relies on the fact that 
most proteins are rigid in comparison to the flexibility of a lipid molecule. Hence, 
when trying to match the hydrophobic region of the protein to the hydrophobic core 
of the bilayer, it is the lipid that will undergo the vast majority of the deforma- 
tion. For the calculations considered here, we assume that leaflet deformations are 
symmetric: whatever happens to the top leaflet is mirrored in the bottom leaflet. 
The deformation is measured as the deviation of the equilibrium position of the 
lipid head-groups by the function u(r) at each position r on the membrane as was 
introduced schematically in Figure [TJ The bending energy takes the form 



G 



(loaf) 



c ) 2 d 2 - 



(14) 



where Kb — 20 knT is the bending modulus of a bilayer (iNiggemann et all . 11995 



Rawicz et all . 120001 ) . equal to approximately twice the bending modulus of a leaflet, 
and the spontaneous curvature of the leaflet, c D , characterizes the leaflet's natural 



tende ncy for a curved state at a hydrophobic-hydrophilic interface (iDan and Safranl . 
19981 ). For many bilayer fo rming lipids , such as phosphatidylcholines, the spon- 
taneous curvature is small (IBoall . 120021 ). In addition to bending, matching the 
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hydrophobic regions of the protein and bilayer necessarily means the bilayer will 
change in thickness, giving rise to a bilayer energy penalty of the form 



2 / V / 



where Z ~ 1.75 nm is the leaflet hydrophobic thickness, and due to membrane volume 
conservation, the bilay er area stretch mod ulus, Ka — 60/«BT/nm 2 , is associated 
with this deformation ( Rawicz et al . 2000l ). Yet another contribution to the free 



energy of deformation in those cases where the membrane thickness is perturbed 
is a local change in the area per lipid as the bilayer thickness varies around the 
protein. Membrane volume conservation arises because the membrane is roughly 



forty times more resistant to vo lume change than area change (ISeemann and Winter! . 



20031 ; iTosh and Collinga . Il986l ) . As a result, if a transmembrane protein locally thins 
the bilayer, lipids will suffer an area expansion in a way that conserves volume. 
Similarly, if the protein locally thickens the bilayer, the area per lipid will locally 
decrease. This implies that the area change near the protein is proportional to the 
compression u(r), and the work done on the bilayer is the integrated area change 
multiplied by tension 

Gt t] = rj^d\ (16) 

where r is the externally applied bilayer tension. Hence, u less than zero corresponds 
to a reduction in the energy of the loading device. All of these contributions can 
be added up to yield the free energy cost associated with thickness variations of the 
two leaflets that can be written as 

In elastic models of this type, the protein dictates the degree of hydrophobic height 
mismatch, u(R) = u , and the angle at which the leaflet contacts the protein at 
the interface between the protein and the surrounding lipids. Far from the protein, 
we expect the bilayer to be flat and slightly thinner in accordance with the applied 
tension, i.e. |Vw(oo)| = and u(oo) = —tI/Ka, respectively. In the case of a 
cylindrical protein we m ake the furth er simplifying assumption that the angle is 



zero (i.e. \Wu(R)\ = 0) (iHuana . Il986l ). The hydrophobic mismatch itself depends 
on membrane properties; changes in membrane thickness are linearly related to the 
hydrophobic mismatch by m = d/2 — I, where d is the hydrophobic thickness of 
the protein. Unlike midplane deformation, the length scale at which the leaflet 
returns to its unperturbed state, A, depends only on fixed elastic parameters of the 
membrane given by 

A = ry^V ^ inm. (18) 

The deformation energy due to thickness variation in the surrounding lipids induced 
by the protein can be written in a simple form when the radius of the protein is 



16 



THE ROLE OF LIPID BILAYER MECHANICS IN MECHANOSENSATION 



larger than A (which is the case for MscL) as 

«(M (B , T)= . Kt (- +7 I_i) 2 ( 1 + V2f). (19) 
The deformation energy scales linearly w ith protein radius and depe nds quadrat- 



ically on the hydrophobic mismatch, u Q (1 Wiggins and Phillipsl . 120051 ). making the 



overall deformation energy particularly sensitive to the hydrophobic mismatch, and 
hence leaflet thickness I □. The deformation energy is fairly insensitive to changes in 
stretch stiffness, Ka (i.e. most terms in the energy are sublinear), and generally in- 
sensitive to changes in the bending modulus since G oc k^ 4 . Additionally, given the 
actual values of the elastic parameters, one finds that the leaflet free energy scales 
roughly linearly with tension, due to the very small value of rj Ka- Like midplane 
deformation, we see that thickness deformation prefers a smaller protein radius. On 
the other hand, in the midplane case, tension always increases the deformation en- 
ergy around a channel while in the case of lipid bilayer thickness variations, the 
tension can either increase or decrease the deformation energy depending on the 
sign of the hydrophobic mismatch. In fact, since the hydrophobic mismatch can be 
either positive or negative {i.e. the protein can be thicker or thinner than the bi- 
layer), tension will increase the deformation energy around a protein that is thicker 
than the membrane {e.g. the closed state of MscL) and decrease the deformation 
energy around a protein that is thinner than the membrane {e.g. the open state), 
as was shown in Figure [7J 

One of the beautiful outcomes of this simple thickness variation elastic theory 
is that the total free energy as a function of protein radius can be written in the 
simple form 

G{R,t) = G (lcaf) (i?,r)- rnF? , (20) 
membrane loading device 

which is reminiscent of classical nucleation theory and results in free energy profiles 
as shown in Figure El At zero tension, the deformation clearly prefers a smaller 
protein radius, limited only by the steric constraints of the protein structure, which 
means that there is a certain minimum radius that the protein can adopt. As 
the tension increases, the quadratic dependence of the driving force on radius will 
eventually overcome the linear dependence of the deformation energy, leading to a 
preference for the open state (corresponding to larger R). We introduce a "hard 
wall" potential at the open radius which provides a severe energy penalty for radii 
larger than the open state radius and argue that this approximation captures the 
idea that opening the channel any further would lead to some degree of energetically 
costly denaturation. This model also captures the correct scale for the critical gating 
tension which is on the order of 1 ksT/nm 2 . 



"''The concept of hydrophobic mismatch is valid when the hydrophobic regions of the protein and 
the bilayer strongly interact, however, this concept has its limits based on the chemistry betwee n 
the lipids and the transmembrane region of the protein {[3, 120031 l2005t iMarkin and Sachd. l2004 . 



and e ventually this condition will be broken if the mismatch is too large (jwiggins and Phillipsl . 
l2005h . 
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Figure 8: Thickness deformation and tension induced energy of a MscL-like channel. 
Competition between the cost of deforming the lipid surrounding a protein and the 
benefit of opening a pore under tension leads to a bistable switch. At zero tension, 
the cost of deformation favors a small protein radius, limited only by the steric 
constraints of the protein. As tension increases, the benefit to opening a pore is 
comparable to the energetic cost to deform the lipids surrounding the protein, and 
a larger protein radius is now possible. At high tension, the potential energy of 
the loading device far outweighs the deformation cost and a larger protein radius is 
favored, again limited by the steric constraints of the protein. 
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It is of interest to compare the energy scale implied by this elastic model to mea- 
sured values. The free energy chan ge of MscL gat i ng was measured to be g : 51 k^T 
using native bacterial membranes (lAnishkin et all . 120051 ; IChiang et all . 120041 ) . If one 
uses the independently measured geometrical properties of the channel, contained 
in Table [fl and elastic properties of pure bilayers (in the text) to calculate the free 
energy of the closed and open states, their difference is approximately 55 ksT at the 
critical tension of 2.5 /csT/nm 2 . Though very encouraging, this close correspondence 
depends sensitively upon the choice of hydrophobic mismatch, as dictated by the 
channel structure and bilayer thickness. 



3.3 Approximating Bilayer Deformation: The Variational 
Approach 

In previous sections, we performed cursory derivations of the energy functionals 
which govern membrane shape for both midplane and membrane thickness defor- 
mations. In order to extract meaning from these energies, we had to minimize the 
free energy functionals of eqns. [11] and [17] with respect to membrane shape. To 
solve the full problem, the conventional scheme (used to obtain the earlier quoted 
results) is to use the calculus of variations to derive a corresponding partial dif- 
ferential equation in the unknown deformation fields h(r) and u(r). A useful and 
intuitive alternative is to adopt a variational approach in which we guess a family of 
solutions (called 'trial functions') that depend upon a small set of parameters and 
then minimize the deformation energy with respect to those parameters. 

For simplicity, we will showcase this method for one- dimensional membranes 
which amounts to the approximation that the protein radius is larger than the 
natural length-scale of deformation, schematized in Figure |9] We will use the varia- 
tional approach to find an approximation for the functions h(r) and u(r) with their 
related energies, and in the process derive the natural length-scale of deformations 
in both cases. Picking a 'good' trial function is intimately related to the success 
of the variational approach. The choice of the trial function is often dictated by 
what we know about the character of the solution. In this case, we know that in 
the near-field the protein is locally disturbing the bilayer by inducing bending or 
hydrophobic mismatch. In the far-field, these disturbances should decay back to a 
flat bilayer. Keeping in mind that most of the energy cost is stored in the local 
disturbance around the protein, we want a trial function that has locally varying 
character around the protein and then a simple decay far from the protein. Such 
a trial function (call it f(x)) could be constructed using a local disturbance, g(x), 
within a decaying envelope 

f(x) = g(x)e-*'\ (21) 

The constant A is an as-yet undetermined natural length scale of deformation and 
will emerge from the minimization process itself. Further, this choice of an exponen- 
tial envelope essentially guarantees that the membrane returns to its unperturbed 
state far from the protein. 

As a practical tool for calculation, our choice of g(x) should have enough param- 
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Figure 9: Protein-induced line tension. Deformation of the membrane around an 
ion channel can be described using a line tension. This line tension is obtained by 
solving for a one-dimensional deformation energy per unit length and then imposing 
that energy around the circumference of the channel. The diagrams above show 
the sequential wrapping of this one-dimensional line tension around a cylindrical 
channel. 

eters to reproduce the given boundary conditions. In addition, we want to choose 
g(x) such that the free energy is a simple function of these parameters. The power 
of the variational approach is that once we have written the energy in terms of these 
variational parameters, the best version of f(x) is, by definition, the one that min- 
imizes the energy. Thus, for instance, if the trial function has two free parameters 
a and b, f(x;a,b), finding the best trial function amounts to solving a system of 
algebraic equations defined by 

— G[f(x; a, b)} = and ^G[f{x; a, b)} = 0, (22) 

where the brackets indicate the energy, G, is calculated using the trial function /. 
This variational strategy can also be used as the basis of numerical approaches in 
which the membrane deformation is represented using finite elements, for example. 
In this case, the trial functions permit us to determine the energy to an arbitrary 
degree of accuracy. Our strategy in the remainder of this section is to use the 
simplicity of the variational approach to find approximate energies for the midplane 
and thickness deformations imposed by membrane proteins. 

3.3.1 Variational Approach for Midplane Deformations 

Our goal is to obtain an approximate expression for the one-dimensional energy due 
to midplane bending given by 

G <-»> = MjT (^M.)) 2 + f (^M.)) 2 ) P3) 

The presence of the 2nR in this expression is due to the fact that we are computing 
the energy per unit length for a deformed bilayer, as shown in Figure [HI and must 
then multiply by the length (the circumference) of deformed material. 

The strategy employed in the variational approach is to plug the trial function 
into the free energy functional and compute the resulting energy, which depends 
upon the parameters in the trial function. Our trial function has the form 

h(x) =g{x)e~(x). (24) 
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The choice of g(x) can be made based upon the boundary conditions. In particular, 
at the boundary of the protein, we require that 

A-h(x)\ x=0 = 9, (25) 

which tells us that we can make the choice g(x) = constant. Applying this boundary 
condition yields the functional form 

h(x) = -0Ae~(*), (26) 

where the only remaining undetermined parameter is the length scale, A. This trial 
function can be plugged into eqn. [23] and the integral is easily evaluated to yield the 
free energy 

G^\\) = ^R6 2 (T\ + ^y (27) 

The next step in the variational strategy is to minimize the free energy with respect 
to A, 

|lG (mid) (A) =0 _> \ = J^, (28) 

which upon substitution yields 

G ( m id) = q\ k & = e 2 nR^¥. (29) 
A 



This is precisely the asymptotic (R^t/kj, > 1) form of eqn. [T2lfor midplane bending 
energy, and our minimization correctly defines the natural length-scale of midplane 
deformation. 



3.3.2 Variational Approach for Membrane Thickness Deformations 

A similar analysis can be made for the one- dimensional deformations induced by 
hydrophobic mismatch. In this case, the free energy functional in the absence of 
tension can be written as 

f (W) ,+ ?0>)> ™ 

We adopt the same functional form for the trial function, namely, 

u(x) = g(x)e~^). (31) 

In this case, we specify two boundary conditions in the near-field; there is a hy- 
drophobic mismatch which demands 

u(R) = u , (32) 
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Figure 10: Comparison of the exact and variational solutions for the thickness de- 
formations around a protein. The variational approach generates an approximation 
to u(x) which is close to the exact solution. The protein radius is R/X = 3. 



and the leaflet has a particular slope at the membrane interface, which we will set 
to zero, 

-^-u(x)\ x=0 = 0. (33) 

In order to accommodate these two boundary conditions, g(x) must have two free 
parameters. As a result, we pick the simplest function which has two degrees of 
freedom, namely a line, and hence set g(x) = ax/X + b, where a and b are constants. 
Applying the two near-field boundary conditions constrains the trial function to the 
form 

u(x)=u (l + ~)e-(*), (34) 

where A is a free parameter with respect to which the energy must be minimized. 
Using this trial function, the free energy can be written as a simple expression of 
the form 

G M(\) = *K b ulR(l^\ + ±y (35) 
Minimizing the energy with respect to A gives 

£g*-»(A) = -> A=nnV#V (36) 



dX w V 5 / 
which upon substitution gives the membrane thickness energy 



G (lcaf) = ( I V {% ) '™ h ulR. (37) 



3/ \Kbl 

Again, the variational approach has reproduced the correct asymptotic form of the 
energy with a small multiplicative error (see eqn. IT5]) ; the exact asymptotic result 

3 

has v^2 instead of (|) 4 , introducing an error of ~ 4%. 

Finally, there are many forms of u(x) which yield roughly the same energy, but 
how does the exact deformation shape compare with our minimized trial function? 
Here too, the variational approach gives a trial function that nearly matches the 
exact result as shown in Figure [TUl 
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Figure 11: Contributions to the free energy. This figure shows how the different 
modes of deformation contribute to the overall free energy budget of the membrane- 
protein system. The energies are written asymptotically to show their dominant 
scaling with the relevant parameters. For the sake of simplicity, we did not address 
how spontaneous curvature factors into the free energy budget. However, a thorough 
discussion of both midplane and leaflet s pontaneous curvature energy contributions 
are found in IWiggins and Phillips! (120051 ) . 



3.4 Distilling the Design Principles 

Having explored how midplane bending, thickness variation and area change are cou- 
pled to tension and the geometric features of the MscL channel, can we distill general 
rules for what makes a membrane protein mechanosensitive? One simple statement 
is that under tension an increase in protein area is always favored, regardless of 
bilayer elastic properties, because an increase in area lowers the potential energy of 
the loading device. Conversely, both midplane and thickness deformations prefer a 
smaller channel, because a larger radius results in a larger annulus of deformed lipid 
and hence a larger free energy penalty (except in the case where the spontaneous 
curvature favors a larger radius (I Wiggins and Phillips! l2005h ). With the area change 



preferring a larger radius and deformation preferring a smaller radius, we have the 
necessary energetic competition that ultimately leads to bistability. This also means 
the sign of the free energy change due to deformation (midplane or thickness) must 
be positive. Hence, the channel is going from a closed state with less deformed lipid 
surrounding it, to an open state with more deformed lipid surrounding it. The con- 
tributions to the free energy budget of a mechanosensitive protein, like MscL, due 
to channel-area change and membrane deformations are shown in Figure [TTJ The 
basic point of this picture is to show how various contributions to the free energy 
scale with the radius (R) and the elastic parameters. 

Midplane deformation is the deformation which depends most simply on mem- 
brane properties, since it is only linked to the bending modulus. Additionally, its 
tension dependence is such that the cost of the deformation always increases with 
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tension and angle, hence we know that tension in addition to preferring a larger 
protein, also wants a more cylindrical protein in the case of midplane bending. This 
allows us, within the limitations of our theory, to put an upper bound on the cost 
of midplane deformations. Taking the lytic tension as an upper bound, a nominal 
bendi ng modulus of 20 / crT, and 9 = 0.6 as a reasonable value of the membrane 
slope ( ITurner and Send . l2004h . the maximum energetic cost of deformations for a 
protein of radius R (in nm) is ~ R x 9 fc^T/nm. 

Thickness deformation depends on all the elastic parameters; bending modulus, 
area stretch modulus, and membrane thickness. The tension dependence of thickness 
deformation energy is also more complex, though a general principle does emerge. 
We know that tension can increase or decrease the overall thickness deformation 
energy, but the general principle is that it always prefers the protein to have the 
same hydrophobic thickness as the bilayer, though the bilayer thickness is itself de- 
creased as tension increases. The other important feature to note is that a decrease 
in the thickness of a transmembrane protein is always accompanied by an increase 
in the area of the membrane surrounding the protein due to volume conservation of 
the membrane. This change in membrane area is indistinguishable from a change 
in protein area. Indeed for MscL, the measured area change is probably a mix of a 
change in the areal footprint of the protein, and a local increase in the membrane 
area surrounding the protein, together giving the measured value of ~ 20 nm 2 . An 
estimate of the upper bound of leaflet deformations is made by assuming the maxi- 
mum u Q = 0.5 nm, then the maximum change in free energy for a protein of radius 
R (in nm) is ~ R x 22/c^T/nm at zero tension with the given elastic parameters 
(see Tabled]). This illustrates that while both midplane and thickness deformations 
are important factors in determining the preferred protein conformation, thickness 
deformations are generally associated with a slightly higher energy scale. 



4 Experimental Considerations 



Much of our knowledge of the function of mechanosensitive channels, including 
MscL, comes from detailed electrophysiology studies where gating of the chan- 
nel is monitored by sharp differences in the ion flux through a membrane patch 



net is monitored by snarp dirterences in tne ion mix tnrougn a membrane patci 
flAnishkin et all . 120051 : IChiang et al 12004 IPerozo et al l2002bl ; ISukharev et all . 12001 



19971 . I1999T ). A small voltage (~ 50m V) is applied across a patch of membrane at the 
tip of a micropipette. As a function of pressure difference, channel opening events 
are recorded as stochastic changes in patch current by an ammeter with picoamp 
(pA) sensitivity. This truly amazing single-molecule spectroscopy technique allows 
the experimenter to adjust the voltage as well as the pressure difference across the 
membrane as shown in Figure [12j The pressure difference across the membrane 
translates into a lateral membrane tension (via the Laplace- Young Relation), re- 
sponsible for gating the mechanosensitive channel. However, there are two serious 
problems with this method when probing the mechanisms of mechanosensitive chan- 
nels. 

Arguably, the most serious problem is that often pressure difference (J/m 3 ) across 
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the membrane is taken to be the input variable of prime importance, when in fact 
tension (J/m 2 ) is the membrane parameter which governs mechanosensitive gat- 
ing. Pressure difference is linearly related to tension via the radius of curvature 
of the membrane, hence in principle the fix is straightfor ward - image the mem- 
brane patch (see Figu re [T2]) . While certainly not impossible ( iMoe and Blountl . 12005 



Sukharev et all . Il999l ). the membrane patch can be difficult to ima ge due to its small 



size a nd the fact that it is inside the micropipette. A recent study (IMoe and Blount 



20051 ) demonstrated the importance of measuring tension in lieu of pressure differ- 
ence. It was shown that using the standard methods for creating "identical" mi- 
cropipettes, the measured characteristics of a channel varied significantly. However, 
when the membrane patch was imaged and tension used as the principle input vari- 
able, the same data collapsed to within a few percent of each other. In general, if 
one could perfectly control the size and shape of the micropipette tip used for con- 
tacting and sealing the membrane patch, all measurements would be related by a 
single constant (the radius of curvature). However, variations in micropipette shape 
and size, as well as variations in how the membrane contacts the pipette tip all lead 
to potentially large variations in the perceived gating characteristics of the channel. 
Additionally, it is difficult to compare the wealth of quantitative data coming from 
electrophysiology studies to theoretical models when pressure difference, instead of 
tension, is used as the principle inp ut variable. Tension is routinely measured in mi- 
cropipette aspiration experiments (IRawicz et all . |2000| ). and in fact, single- c hannel 
electrophysiology recordings are possible in such a setup (IGoulian et all . Il998l ) using 
ion channels with conductances lower than MscL. Hence, this technique might pro- 
vide a useful way to apply known membrane tension to reconstituted MscL channels 
in well characterized membranes. 

With tension being used as the variable of prime importance, electrophysiology is 
poised to put the continuum mechanical view to the test, elucidating the role of lipids 
in ion chan nel function. In par ticular, the elastic properties of many lipids have been 
measured ( Rawicz et al . 2000l ). enabling a careful examination of the dependence of 
gating energy on lipid carbon chain length. The simple continuum view we set forth 
here predicts a quadratic dependence of the lipid thickness deformation energy on 
hydrophobic mismatch, which is directly linked to carbon chain length. This, of 
course, has implications for both the function of various transmembrane proteins, 
and comments meaningfully on the ability of bilayer thickness to segregate proteins 
in biological membranes. 

A second class of intriguing experiments concerns the mechanosensitivity of 
other ion channels and recepto r s, gen erally regarded not to be mechanosensitive 
(jCalabrese et all . 120021 ; iGu et all . 1200 ll ). This is both interesting from a functional 
standpoint, in an effort to understand the full physiological effects of these pro- 
teins, and as a tool for understanding structural features such as the motions of 
transmembrane helices. Performing a similar experiment where lipid carbon chain 
length is varied around a voltage-gated ion channel (for example) could reveal hidden 
mechanosensitivity, and energetic analysis from such an experiment could comment 
on the degree of height and area change during the gating transition. 

The second problem facing a complete understanding of the function of 
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mechanosensitive channels is that for many such channels volumetric flow, and not 
ion flux, is the relevant physiological parameteiQ. Hence, ion flux is used as a surro- 
gate measurable in place of the true physiological output of the channel. One could 
argue that ion flux is proportional to volumetric flow, however this assumes that the 
way ions flow through the channel pore is identical to the way water flows through the 
pore. Ex periments have elucidated the roughly ohmic nature of mechanosensitive 
channels ( ICruickshank et all . 119971 : iPerozo et a ], l2002bl ) at low voltage (< 80m V), 
however we know essentially nothing about how a pressure gradient across the mem- 
brane translates into a volumetric flow. Even the simplest continuum approximation 
(Hagen-Poiseuille flow) would predict a non-linear function relating the area of the 
channel pore to the volumetric flow, in contrast to the (roughly) lin ear re l ation ship 
between ion flux and channel pore area as predicted by Ohm's Law ( IHillei . Il968l ). It 
would be of considerable physical and physiological interest to expand our under- 
standing of fluid flow at the molecular level, by measuring the relationship between 
pressure gradient and volumetric flow through a large-pore channel like MscL. 



5 Cooperativity and Interaction between Trans- 
membrane Proteins 

One intriguing consequence of the deformations induced in membranes by ion chan- 
nels is that channels will interact. These interactions can lead to cooperativity in the 
gating of neighboring channels and can also induce spatial ordering of the proteins. 
These interactions can be thought of as arising from two different effects: those of 
elastic origin and those of thermal origin. The elastic forces are purely an enthalpic 
effect coming from a minimization of the deformation energy around two proteins 
separated by a given distance. The thermal forces are entropically driven by the 
thermal fluctuations of the membrane and are analogous to Van der Waals forces. 



5.1 Enthalpic Interactions 



As discussed above, proteins which change the membrane thickness or bend the 
membrane midplane produce deformations which extend anywhere from a few 
nanometers (thickness) up to tens of nanometers (midplane) from the protein edge. 
As two proteins approach each other, their respective deformation fields overlap re- 
sulting in a deformation profile between them that is different than either of them 
produce separately. In this case, the total deformation energy of the system is de- 
pendent on the separation between the two proteins and results in an interaction 
potential which is dependent upon the conformation of the proteins. These forces 
arise purely from the mechanical attributes of the deformed membrane and have no 



I The mechano sensitive bacterial channel MscS (jBass et al . 20021 Levina et al . 1999; 



Pivetti et all . l2003h is another example. Although, there are also mechanosensitive channels that ap- 



pear to be highly ion selective, such as the bacterial mechanosensitive i on channel MscK dLi et all . 
2002f) and the K2P family of mammalian mech anosensitive channels ( Franks and Honorei 20041 : 
Lauritzen et~al l2005t iMaingret et al 1 1 999l [2000h . 
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Figure 12: Measurement of tension vs. pressure difference in an electrophysiological 
experiment. A channel protein (small blue rectangle) is embedded in a membrane 
patch (green). A potential of roughly 50mV is applied across the sealed membrane 
patch, and channel opening events are measured by an ammeter (circle) with pi- 
coamp (pA) sensitivity, a) At low pressure difference, the tension in the patch is 
low, the mechanosensitive channel is in the closed conformation, and the patch has 
a very large radius of curvature. The plot to the right shows normalized channel 
current as a function of time for a simulated channel; the open state has low oc- 
cupation at low tension, b) At high pressure difference, the tension in the patch 
is high, the mechanosensitive channel will occupy the open state, and the radius of 
curvature (r) is on the order of microns. The plot to the right shows the open state 
has high occupation at high tension, c) Optical micrograph of vertically oriented 
membrane patch at low (top) and high (middle and bottom) pressure differences, 
illustr ating the decreas e in th e radius of curvature with increase pressure difference 
(from iMoe and Blountl fl2005h K The scale bar is 5/im. 
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entropic component. We know that midplane and thickness deformations are inde- 
pendent, and hence there are distinct interactions due to midplane and thickness 
deformations. 

Pairwise interactions due to midplane deformation using eqn. [11] have previ- 
ously been calculated for a variety of membra ne curvature environments and pro- 
tein shapes at zero tension (IChou et all . 120011 ). Using a bilayer bending modulus 
of ~ 100 ksT, attractive interactions of order ~ 1 — 5 /c^T were found when the 
proteins were separated by 1 - 2 protein radii (which we estimate to be 5 - 10 nm 
measured center-to-center for a typical transmembrane protein). If we adjust the 
energy scale to be consistent with a phosphatidylcholine bilayer bending modulus of 
~ 20 kgT this lowers the interaction energetics to ~ 0. 5 — 3 k^T. These int eractions 
tend to be long-ranged with a power- law decay of 1/r 4 ( IGoulian et all . Il993l ) . Simple 
pairwise interaction will be inadequate to describe the nature of interactions between 
more than two proteins. This arises because one protein can shield other proteins 
from feeling the deformation of a neighboring protein, and hence interactions are 
not (in general) pairwise additive. Apart from direct numerical simulation, there 
are few analytical (theoretical) tools which allo w one to study how m any interacting 
proteins in close proximity behave as a group (IHarroun et all . Il999l ). 



Like midplane deformations, the thickness deformation fields extending from 
the edges of two proteins will overlap and interact as the proteins come into close 
proximity (lAranda-Espinoza et all . Il996l : iDan et all . Il993l ). We provided evidence 
that lipids likely influence the function of MscL through thickness deformations 
and once again we will appeal to MscL study for interacting membrane 

proteins. The short-range nature of thickness deformations (essentially exponential 
decay) means there is no power-law asymptotic formula for their interaction, though 
we numerically explored these interactions for all possible conformations of two MscL 
proteins as shown in Figure [T31 As we saw with single proteins, the energetic scale 
of thickness interactions is generally higher than with midplane deformations, and 
can vary greatly depending on the hydrophobic mismatch. The leaflet interactions 
between two MscL proteins are appreciable when they are within several nanometers 
of each other, and ranged from ~ 2 — 25 k%T depending on the protein conformations 
and the tension in the membrane. This kind of s hort-ranged interaction might play 
an important role in mem brane protein function (jBotelho et all . 120061 ; iGoforth et all . 
20031 : Molina et all . l2006f ). given the nominal density of transme mbrane proteins i n 
biological membranes leads to spacings on the order of 10-100 nm (IMitra et all . 12004 ). 

Additionally, the interactions due to thickness variations can be either attractive 
or repulsive depending on the shape of the proteins. The general principle that 
emerges is that 'like' proteins attract and 'unlike' proteins repel (in contrast to 
electrostatics) as shown in Figure [131 Proteins whose hydrophobic mismatch has 
the same sign (i.e. both are taller or both are shorter than the membrane) lead to 
net attractive interactions; proteins with opposite signs of hydrophobic mismatch 
lead to repulsive interactions. Later in the article, we will demonstrate that these 
conformation-dependent interactions can communicate state information between 
two proteins, leading to cooperative channel gating. 
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Figure 13: Conformation-dependent interactions between two MscL channels. As 
two MscL channels (blue) come close to each other, regions of deformed lipids (green) 
overlap leading to deformation induced interactions, a) Deformation surrounding 
two closed MscL channels, b) Deformation surrounding a closed and an open MscL 
channel, c) Deformation surrounding two open MscL channels. The relative sizes of 
the open channel, the closed channel, and the lipids are roughly correct, d) Inter- 
action potentials for the three configurations shown in a, b and c. External tension 
weakens the interaction between two open channels (V ao ) and strengthens the inter- 
action between two closed channels (V cc ), but has almost no effect on the interaction 
between an open and closed channel (V oc ). The open-open and closed-closed interac- 
tions are both more strongly attractive than the open-closed interaction, indicating 
that elastic potentials favor interactions between channels in the same state. The 
'hard core' distance is where the proteins' edges are in contact. 



5.2 Entropic Interactions 



A second class of forces between membrane proteins arise due to membrane fluctu- 
ations. Like most entropic forces, the thermal interactions between transmembrane 
proteins are fairly weak, on the order of a few ksT. Two fluctuation-induced forces 
have been studied in some detail in the lite rature; a long-ranged Casimir force due to 
the s urface fluctuations of the membrane (IGoulian et all . Il993t iPark and Lubenskyl . 
1990) , and a very short-ranged depletion fo rce due to the excluded volume of lipid 
molecules between two membrane proteins (ISintes and Baumgartnerl . Il997t ) . 

The Casimir force between two membrane proteins arises because the available 
spectrum of fluctuations of the membrane-midplane depend on the distance be- 
tween two proteins. Entropically, the membrane-protein system seeks to maximize 
the number of available modes of fluctuation and hence an energetic potential exists 
between two transmembrane proteins in a fluctuating, thermally active membrane. 
Through a series of elegant calculations, this force was shown to have a 1/r 4 asymp- 
totic form, wh e re r i s the center-to-center distan ce between two cylindrical proteins 
( Goulian et ah 1993 ; Park and Lubenskyl . 19961 ). If we presume that it is approxi- 
mately correct for small separations (r ~ 2i?), this implies an attractive potential 
with an energy scale of ~ 1 k B T. 
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Figure 14: Conformational statistics of interacting MscL proteins. Interactions be- 
tween neighboring channels lead to shifts in the probability that a channel will be 
in the open state (dashed lines). The sensitivity and range of response to tension, 
dP opcn /dr, are also affected by bilayer deformations (solid lines). P opan and dP open /dr 
are shown for separations of 0.5 nm (red) and 1.5 nm (green) with reference to non- 
interacting channels at d = oo (blue). Interactions shift the critical gating tension 
for the closest separation by ~ 12%. Additionally, the peak sensitivity is increased 
by ~ 90% from ~ 5nm 2 /ksT to ~ 9.5nm 2 //cbT. 



Lateral density fluctuations of lipids in the membrane also lead to entropic forces 
between proteins. Using Monte Carlo simulations, these entropic depletion forces 
(also called 'excluded volume forces') between cylindrical proteins were shown to be 
appreciable only when the prote ins' edges were within ~ 1 lipid molecular diameter 
(jSintes and Baumgartnerl . 119971 ). For cylindrical proteins, with diameters on the 
order of ~ 1 — 2 nm, direct edge contact resulted in a favorable interaction with an 
energy scale of ~ 2 ksT. 



5.3 Protein Conformations Affected by Interaction 

As noted above, the elastic interactions between ion channels such as MscL depend 
upon protein conformation. In earlier sections of the paper, we established that the 
equilibrium conformations of the channel are entirely determined by the free energy 
difference between the two states. As a result, elastic interactions which change the 
energy of a two-channel system will affect the probability that we measure any one 
channel in the open state. In fact, electrophysiology (see Figure [£2]) is well suited to 
such measurements where the total amount of time spent in the open state divided 
by the total measurement time is the open probability. 

The free energy difference between th e open and closed states of a single channel 



is roughly 50 ksT (IChiang et all . |2004| ). which implies the energy scale for two 
channels is roughly 100 k^I '. We have also seen that two MscL channels in prox- 
imity have interactions on an energy scale of roughly ~ 20 ksT as shown in Figure 
Two open channels have a strong, favorable interaction that can significantly 
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alter the open probability of a given channel relative to the isolated channel value 
as shown in Figure dH Such interactions also affect channel 'sensitivity', defined 
by the derivative of the P opon curve with respect to tension, which quantifies how 
responsive the channel is to changes in the driving force, in this case tension. The 
full-width at half maximum of this peaked function is a measure of the range of 
tension over which the channel has an appreciable response. In general, the area 
under the sensitivity curve is equal to one, hence increases in sensitivity are always 
accompanied by decreases in range of response, as demonstrated by the effects 
of the beneficial open-open interaction on channel statistics. The critical gating 
tension and sensitivity are essentially the key properties which define the transition 
to the open state, and are analogs to the properties which define the transition 
of any two-state ion channel. Hence, the elastic interactions can affect channel 
function on a fundamental level. 

Conclusion 

The goal of this article is to take stock of the role of lipid bilayer deformations in 
mechanosensation. More precisely, we have argued that the lipid bilayer is not a 
passive bystander in the energetics of channel gating. As a result, by tuning mem- 
brane properties it is possible to alter channel function. We have emphasized two 
broad classes of membrane deformation that are induced by the presence of a trans- 
membrane protein: i) deformation and bending of the midplane of the lipid bilayer, 
ii) variations in the thickness of the lipid bilayer that are induced by hydrophobic 
mismatch. As a result of these deformations, there is a free energy cost to changing 
the radius of a channel since the open state implies a larger annulus of deformed 
material and hence a higher free energy. This deformation energy competes with 
the energetic relaxation of the loading device. 

One of the key reasons for performing theoretical analyses like those described 
here is that they permit us to sharpen the questions that can be asked about a given 
biological problem. This sharpness is ultimately most meaningful if it is translated 
into precise experimental predictions. The theoretical results described here suggest 
a variety of experimental predictions. 

• Dependence of gating tension on hydrophobic mismatch. Previous work has 
already shown that lipid bilayer tail lengths can alter channel gating by chang- 
ing the hydrophobic mismatch. To more precisely examine this relationship, 
careful measurements of the membrane tension need to be made, as opposed 
to pipette pressures, to elucidate the energetics underlying gating. Alterna- 
tively, mutagenesis could be used to explore the same effect by changing the 
hydrophobic thickness of the protein. 

• Hidden mechano sensitivity in other classes of channels and receptors. The re- 
sults described here have been applied to the case study of MscL. However, we 
argue that any transmembrane protein that varies its radius or hydrophobic 
thickness upon conformational change will exhibit mechanosensitivity. Fur- 
thermore, the way tension affects the function of these proteins might help 
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elucidate the classes of structural changes that occur during their conforma- 
tional change. 

• Cooperative gating of channels. As a result of the elastic deformations induced 
in the lipid bilayer by mechanosensitive channels, nearby channels can commu- 
nicate their conformational state, resulting in cooperative gating. This coop- 
erativity should be observable in electrophysiology experiments as a change in 
the critical tension and channel sensitivity with an increase in channel density. 



Shortcomings of the Theory. Obviously, the use of simple ideas from elasticity 
theory to capture the complex process of mechanosensation provides a caricature of 
the real process. One signature of the shortcomings of this kind of approach is the 
fact that singl e amino acid substitut i ons c an completely alter the properties of cer- 
tain proteins (jYoshimura et all . Il999l . 12004 ) . This serves as a warning of the pitfalls 
of models that ignore atomic-level details and their impact on biological function. A 
second class of complaint that can be registered against the models described here 
is that we have ignored material heterogeneity. In particular, biological membranes 
are built up of a broad range of different lipids and are riddled with membrane 
proteins. As a result, it is not clear if an elastic description like that described here 
is appropriate, and if it is, how to select the relevant material parameters. 

Regardless of the difficulties highlighted above, it is clear that the emergence of 
an increasing number of structures of ion channels coupled with functional studies of 
these proteins has raised the bar for what should be expected of theoretical models of 
channel function. The central thesis of the work described here is that the presence 
of the lipid bilayer provides another way in which these systems can be manipulated. 
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